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ABSTRACT 

We present two-dimensional (2D) stellar and gaseous kinematics of the inner 
'^130x180 pc^ of the Narrow Line Seyfert 1 galaxy NGC4051 at a sampling of 
^ I 4.5 pc, from near-infrared X-band spectroscopic observations obtained with the Gem- 

ini's Near-infrared Integral Field Spectrograph (NIFS) operating with the ALTAIR 
adaptive optics module. We have used the CO absorption bandheads around 2.3 /^m 
to obtain the stellar kinematics which show the turnover of the rotation curve at only 
~55pc from the nucleus, revealing a highly concentrated gravitational potential. The 
QQ ' stellar velocity dispersion of the bulge is « 60 km s~^ - implying on a nuclear black hole 

, mass of K, 10® Mq - within which patches of lower velocity dispersion suggest the pres- 

ence of regions of more recent star formation. From measurements of the emission-line 
profiles we have constructed two-dimensional maps for the flux distributitions, line ra- 
tios, radial velocities and gas velocity dispersions for the H2, Hi and [CaVlii] emitting 
gas. Each emission line samples a distinct kinematics. The Br7 emission-line shows 
no rotation as well as no blueshifts or redshifts in excess of 30kms~^, and is thus not 
restricted to the galaxy plane. The [Caviii] coronal region is compact but resolved, 
extending over the inner 75 pc. It shows the highest blueshifts - of up to — 250kms~^, 
k> \ and the highest velocity dispersions, interpreted as due to outflows from the active 

5_j ■ nucleus, supporting an origin close to the nucleus. Subtraction of the stellar velocity 

' field from the gaseous velocity field has allowed us to isolate non-circular motions 

observed in the H2 emitting gas. The most conspicuous kinematic structures are two 
nuclear spiral arms - one observed in blueshift in the far side of the galaxy (to the 
NE), and the other observed in redshift in the near side of the galaxy (to the SW). 
We interpret these structures as inflows towards the nucleus, a result similar to those 
of previous studies in which we have found streaming motions along nuclear spirals in 
ionized gas using optical IFU observations. We have calculated the mass inflow rate 
along the nuclear spiral arms, obtaining ~ 4 x 10^^ yr^^, a value ^ 100 times 
smaller than the accretion rate necessary to power the active nucleus. This can be un- 
derstood as due to the fact that we are only seeing the hot "skin" (the H2 emitting 
gas) of the total mass inflow rate, which is probably dominated by cold molecular gas. 
From the II2 emission-line ratios we conclude that X-ray heating can account for the 
observed emission, but the II2 A2.1218 /im/Br7 line ratio suggests some contribution 
from shocks in localized regions close to the compact radio jet. 

Key vifords: galaxies: Seyfert - infrared: galaxies - galaxies: NGC4051 (individual) 
- galaxies: kinematics 
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1 INTRODUCTION 

The presence of supermassive black holes (SMBHs) at 
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nowadays widely acc e pted by the astronomy comunity 
l|Gebhardt et al.fl2000l : [Ferrarese fc Merrittll2000l ). Accord- 
ing to this scenario the energy emitted by an Active 
Galactic Nucleus (AGN) is due to the accretion of ma- 
terial onto the SMBH and implies the presenc e of a gas 
reservoir close to the AGN. ISimoes Lopes et alj l)2007h us- 
ing archival Hubble Space Telescope (HST) optical im- 
ages for a large sample of early-type galaxies with and 
without AGNs, found that all AGN hosts have circumnu- 
clear gas and dust, while this is observed in only 26% of 
a pair-matched sample of inactive galaxies. A gas reser- 
voir close to the AGN is also supported by the presence 
of recent star formation in the circumnuclear region of ac- 
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SMBH is occurring include the observation of streaming mo- 
tions in ionized gas along nuclear spirals towards the nucleus 
of nearby active galaxies using two-dimensional (2D) opti- 
cal sp ectroscopy (jFathi et al.ll2006l : IStorchi-BerEmann et al.l 
120071 ). 

Streaming motions as feeding signatures to the nu- 
clear re gion have been previously observed in radio wave- 
lengths. lAdleT^^^^st^fehJ (|l996h . for example, have found 
streaming motions towards the center along the spiral arms 
of M 81 in H I, while Mundell fc Shone (,1999 ) have found 
similar streaming motions towards the nucleus along the 
weak bar of NGC4151. Closer to the center, most of the 
gas is in the molecular phase, and CO observations have 
been used to map the g a s kinematic s and inflows (e.g. 
Garcfa-Burillo et all l2003l : iKrips et all l2005l: iBoone et al] 



20071 ). 



Molecular hydrogen emission lines are also relatively 
strong in the near-IR K-h&nd spectra of active galax- 
ies, and previous studies suggest that its distribution 
and kinematics is distinct from that observed in the 
other emiss ion lines, which are usually d ominated by out - 
flows (e.g. ICrenshaw fc Kraemerl I2OOOI : iDas et ah I I2OO5I ). 
In the Seyfert galaxies NGC 2110 a nd Circinus, for ex- 
ample, IStorchi-BerEmann et al.l l|l999l) have found broader 
emission-line profiles for [Fell] and Pa/3 than for the 
H2A2.1218/im emission-line in a study using near-IR long- 
slit observations, suggesting a stronger influence of nuclear 
outflows on the former emission lines and a different origin 
for the H2 emitting gas, consistent with the colder kinemat- 
ics of the galaxy disk. More recently, using two-dimensional 
(hereafter 2D) nea r-IR spectroscopy of the Seyfert galaxy 
ES0 428-G14, we (jRiffel et al.l |2006l) have found that the 
H2 emission distribution was mostly restricted to the plane 
of the galaxy and was less affected by the AGN outffow than 
the [Fell] and Pa/3 emission lines. 

In this paper we use adaptive optics IFU spectroscopic 
data obtained with th e Gemini's Near-infrare d Integral Field 
Spectrograph (NIFS - iMcGresor et al.ll2003l '). in the near-IR 
K-hanA at a sampling of O'.'lxO'.'l, to investigate the stellar 
and gaseous kinematics of the inner ~ 3 x 4 arcsec^ of the 
nearby active galaxy NGC 4051. The stellar kinematics will 
be used to constrain the galaxy potential and mass of the 
SMBH, but our main goal is to look for signatures of feeding 
mechanisms at parsec scales through the H2 kinematics. 

NGC 4051 is a SABbc galaxy at a distance of only 



■-^9.3 Mpc (|Barbosa et al.]|2006l ). such that I'.'O corresponds 
to 45 pc at the galaxy. It harbors one of the closest AGN, 
classifled as a Narrow-Line Seyfert 1 (NLSl). We have se- 
lected NGC 40 51 for this study part ially on the basis of the 
recent work of iBarbosa et al.) (|2006| ). who obtained 2D stel- 
lar kinematics using the IFU of the Gemini Multi-Object 
Spectrograph (GMOS) to observe the stellar absorption lines 
of the Calcium triplet around 8500 A. These authors have 
shown that the turnover of the stellar rotation curve occurs 
at only R ~50 pc from the nucleus, indicating that the stel- 
lar motions are dominated by a highly concentrated gravita- 
tional potential. As NIFS with the ALTAIR adaptive optics 
module provides a much better image sampling and resolu- 
tion than the GMOS-IFU, we decided to further investigate 
the kinematics of the nuclear region of NGC 4051 in order 
to better constrain the gravitational potential and the mass 
of th e SMBH. In adition, this galaxy shows stro ng H2 emis- 
sion (|Riffel. Rodriguez- Ardila fc Pastorizall2006l ). making it 
a good candidate for a study of its kinematics. 

Previous studies of NGC 4051 include HST narrow- 
band [O III] images which show an unresolved nuclear 
source and faint extend ed emission (by l'.'2) alon g the po- 
sition angle PA=100° (|Schmitt fc KinnevI Il995 ). the ap- 
proximate direction of the align ment of two radio compo - 
nents a1. 6 cm separated by 0'.'4 (|Ulvestad fc WilsonI [l983 ) . 
IVeilleux] (|l99ll ) reported that the proflles of optical forbid- 
den emission lines present blue wings reaching velocities of 
up to 800kms~^, and proposed a model for the narrow- 
line region (hereafter NLR) with outflows and obscuring 
dust. Evidence for o utflow s have also been observed by 
IChristopoulou et al] (Il997l) for the [Qui] emission to l"5 



NE of the nucleus. iNagao. TaniEUchi fc Muravamal (|2000l ) 
using the [FeX]A6374 emission line report a high ionization 
region extending to 3'.'0 SE of the nucleus. [Lawrence et al.l 
(1985) found a X-ray variability on time scales of the or- 
der of 1 hour and lSalvati et all (|l993l ) reported flux changes 
by a factor of 2 in 6 months observed at 2.2 fim and sug- 
gest that emission by dust reprocessing of the UV radiation 
from the nucleus is an acceptable explanation. iPonti et al.l 
(2006) modelled the nuclear emission by two power law com- 
ponents, one due to the AGN and other due to reflection by 
the accretion disk. 

This paper is organized as follows: in section 2 we de- 
scribe the observations and data reduction. In section 3 we 
present the results for the stellar kinematics. In section 4 we 
present the emission-line flux distributions and in section 5 
we present the results for the gas kinematics. In section 6 
we discuss the results and in section 7 we present the con- 
clusions of this work. 



2 OBSERVATIONS AND DATA REDUCTION 

The IFU spectroscopic data were obtained with NIFS oper- 
ating with the ALTAIR adaptive optics module on the 8-m 
Gemini North telescope in January 2006 under the instru- 
ment science veriflcation program GN-2006A-SV-123. The 
IFU has a square held of view of ~ 3'.'0 x 3'.'0, divided into 
29 slices with an angular sampling of 0'.'lx0'.'04. The ob- 
serving procedures followed the standard Object-Sky-Sky- 
Object dither sequence, with off-source sky positions since 
the target is extended, and individual exposure times of 750 s 
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Figu re 1. Top left: WHT K-band large scale image of NGC 4051 
from iKnapen et al.l | |2003|) . The image has been rotated to the 
same orientation of the NIFS observations. The box represents 
the IFU field of view. Top right: 2.12 fim continuum map from 
IFU spectroscopy. Bottom: Spectra at the positons N, A and B 
marked at the top-right panel with the emission lines and CO 
absorption bandheads identified. 



centered at A = 2.2499 ^m. Two set of observations, each 
with three individual exposures, were obtained at different 
spatial positions; the first one centered at the position 0'.'4 
from the nucleus along the PA= —74° and the second at 
a position offset (y.'5 along PA= 106°. The longest extent 
covered by the IFU observations was oriented along the po- 
sition angle PA= 106° , which corresponds to the orientation 
of the line of nodes derived by modelling of the stellar veloc- 
ity field with a rotating disk ( ,Barbosa et al.i,200 6'V We have 
used the K_G5605 grating and the fiher HK_G0603, which 
resulted in an arc lamp line full width at half maximuum 
(FWHM) of 3.2 A. 

The data reduction was accomplished using tasks con- 
tained in the nifs package which is part of Gemini iraf pack- 
age as well as generic iraf tasks. The reduction procedure 
included trimming of the images, flat-fielding, sky subtrac- 
tion, wavelength and s-distortion calibrations. We have also 
removed the telluric bands and fiux calibrated the frames 
by interpolating a black body function to the spectrum of 
the telluric standard star. The final IFU data cube contains 
1160 spectra, each spectrum corresponding to an angular 
coverage of O'.'lxO'.'l, which translates into 4.5x4.5pc^ at 
the galaxy and covering the spectral region from 1.991 
to 2.425 ^im. The total observed field of view 2'.'9x4'.'0 (ob- 
tained by mosaicing the two set of observations) thus cor- 
responds to a region of projected dimensions 130pcxl80pc 
at the galaxy. 

In the top-left p anel of Fig.[l] we pre sent a large scale 
K-hand image from iKnapen et al.l (|2003l ) obtained at the 
William Hershel Telescope (WHT), The central rectangle 
shows the IFU field of view. The large scale image was ro- 



tated to the same orientation of the IFU observations. In the 
top-right panel we present an image obtained from the NIFS 
data cube for the 2.12 /im continuum emission, obtained 
by interpolation of the continuum under the H2 A2.1218 
emission line. In the bottom panels we present three chara- 
teristic IFU spectra: the nuclear spectrum (position N in the 
continuum map), a spectrum from a location at l'.'2 E of the 
nucleus (position A) and another from 0'.'75 W of the nu- 
cleus (position B). Both spectra correspond to an aperture 
of 0'.'3x0'.'3. The emission lines are identified in the nuclear 
spectrum: the H2 lines at A =2.0338, 2.1218 and 2.2235 ^im, 
the Hi Br7 at 2.1661 ^m and the [Cavill] coronal emission 
line at 2.3211 ^m. The ^^CO and "CO absorption band- 
heads used to obtain the stellar kinematics are identified in 
the spectrum from position A. 



3 LINE OF SIGHT VELOCITY 
DISTRIBUTION 

In order to obtain the line-of-sight velocity distribution 
(LOSVD) we have used the penaliz e d Pix el-Fitting (pPXF) 
method of ICappellari fc EmsellemI Hooi) to fit the stellar 
absorptions present in the if-band spectra. The algorithm 
finds the best fit to a galaxy spectrum by convolving a tem- 
plate stellar spectrum with the corresponding LOSVD. This 
procedure gives as output the radial velocity, velocity disper- 
sion and higher-order Gauss-Hermite moments. The pPXF 
method allows the use of several template stellar spectra 
and to vary the weights of the contribution of the differ- 
ent templates to obtain the best fit, minimizing the tem- 
plate mismatch problem. However, the use of pPXF re- 
quires template s whic h match closely the galaxy spectrum. 
lEmsellem et al.l |2003) present a extensive discussion about 
the use of th e pPXF method and the templates mismatch 
problem and ISilge fc Gebhardd (HoO^) present a discussion 
about template mismatch for kinematic fitting using the CO 
absorption bandheads in the near-IR. 

A high signal-to-noise ratio is required for reliable stel- 
lar kinematic measurements using the pPXF method. In 
previous works the data are usually binned to give S/N ra- 
tio between 40 and 60 over the whole field-of-view. In our 
data, ratios smaller than these are observed only very close 
to the borders of the field (S/N ~ 35 measured blueward 
of the first CO band-head). Our typical S/N ratios are 80 
with maximuum values reaching 120 close to the nucleus, 
thus allowing reliable kinematic measurements without spa- 
tial binning of the data. 



3.1 The stellar templates 

In this study we have selected as template spectra, those 
of the spectroscopic library of late spectral type stars 
observed with the Gemini Near Infrared Spectrograph 
GNIRS) IFU using the grating 111 1/ni m in the K-hand 
Winge. Riffel fc Storchi-Elergmannl [20o3). This library is 
composed of spectra of 29 objects, which include dwarf, gi- 
ant and sub-giant stars with spectral types from F7lll to 
M3lll, observed in the spectral range from 2.24 to 2.42 /xm. 
Of these, 23 stars were also observed on a second set- 
ting extending the wavelength coverage down to 2.15 /^m. 
The spectra of the templates have signal-to-noise ratios 
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larger than 50 blueward of the CO 2-0 first-overtone band- 
head. Detailed discussion about th i s hbra ry is presented by 
IWinge. Riffel fc Storchi-Bergmaniil (j2007h . The spectral res- 
olution of the stellar spectra is 3.3 A - very close to the 
spectral resolution of the NIFS data, thus the library can be 
used to fit the stellar kinematics without any resolution cor- 
rection. We have verified this in some tests using available 
NIFS spectra of a few stars. We opted to use the GNIRS 
library because there are too few available template spectra 
with NIFS. We have chosen 16 stars from the library which 
better match the stellar features of NGC4051. 

In order to investigate the influence of the stellar tem- 
plates on the velocity dispersions obtained we have fitted 
the stellar kinematics using individual stars as templates. 
We observed that the large scale structures in o maps are 
similar for all stars, but the mean a values vary significantly 
- higher EWs in the templates result in lower a values for 
the galaxy. This result evidences the importance of using a 
large library of stellar templates in order to obtain reliable 
velocity dispersion measurements. 



3.2 The effect of the [Cavill] coronal Hne 

The nuclear spectrum of NGC4051 (Fig.[T} sh ows the coro- 
nal em ission line of [Cavill] at 2.3211 /im. iDavies et al.l 
l|2006t ) have shown that this line affects the kinematic mea- 
surements obtained from the CO bandheads. In order to in- 
vestigate the infiuence of this line in our measurements, we 
have chosen two spectral regions to fit the stellar kinematics: 
the first one (2.258-2.372 ^.m) includes the [Cavill] line and 
the second (2.258-2.314 /im and 2.346-2.372 ^m) excludes 
this line (together with the "contaminated" CO bandhead) . 
We observed that the radial velocities and velocity disper- 
sions derived using the first spectral range are higher in the 
nuclear region, where the [Cavill] emission line is present, 
than those obtained by using the second spectral region. We 
thus decided to exclude the contaminated ^^C0 3-l band- 
head from the stellar kinematics fitting in the region where 
the [Cavill] emission line is present (the central O'.'S radius). 
In regions away from the nucleus we used the first spectral 
range which includes the ^^C0 3-l bandhead. 



3.3 The stellar kinematics 

In Fig. [2] we show the fits of the stellar templates to the 
galaxy spectra with the program pPXF for four different 
positions: the nucleus; l'.'2 SE and l'.'2NW of the nucleus, 
where we observe the turnovers of the rotation curve and a 
location at l'.'5 E of the nucleus, almost at the border of the 
NIFS field. We observe that the stellar templates fit very well 
the galaxy spectra at most positions, including regions near 
the border of the IFU field, where the signal-to-noise ratio 
is smaller. The results for the nucleus should nevertheless 
be considered with caution, as the fitting may have been 
affected by dust emission and emission lines present in the 
galaxy spectrum but absent in the stellar templates. 

In Fig. Owe present the resulting stellar kinematics. The 
black regions in this figure are masker regions where the 
signal to noise ratio in the spectra was too low to provide 
a reliable fit. In the top-left panel we show the stellar ve- 
locity field, which shows a velocity range from ~ —40 to 




R (arcsec) 



Figure 4. One dimensional stellar rotation curve for a pseudo- 
slit oriented along the major axis of the galaxy. The points are 
the observed radial velocities and the full line is the modeled 
velocities. 



~ 40kms ^ with the turnover occurring at ~55pc from the 
nu cleus in good agreem ent with the velocity field obtained 
bv lBarbosa et all (|2006l ) from optical IFU data. The mean 
uncertainties in the velocities are ~ lOkms"^. The turnover 
and amplitude of the rotation curve can be more easily ob- 
served in the one dimensional cut of the stellar velocity field 
shown in Fig.|4l The one dimensional rotation curve was 
obtained by the averare of the velocities within a pseudo- 
slit with 1 arcsec width oriented along the major axis of the 
galaxy. 

The stellar velocity dispersion map is shown in the top- 
right panel of Fig.O The a map presents values ranging from 
~35 to ~90kms~^, with mean uncertainties of ~ 8kms~^. 
The bottom panels show the higher order Gauss-Hermite 
moments /13 (left) and /14 (right). These moments measure 
devations of the line profile from a Gaussian: the parameter 
/13 measures asym metric deviations and the /t4 m easure sym- 
metric deviations l|van der Marel fc Franxlll993l ) . The values 
/i3 and /14 vary from -0.15 to 0.15 with mean uncertainties 
of 0.03. The highest values of /13 are observed to SE of the 
nucleus and the lowest values to NW of the nucleus. The 
values are nearly zero over most of the IFU field. 



3.4 Kinematic Modelling 

Since the stellar velocity field is dominated by rotation, it 
was fitted with a velocity model produced by a Plummer 
potential in order to obtain the systemic velocity, orientation 
of the line of nodes, bulge mass and the position of the 
kinematical center. The Plummer potential is given by: 



GM 



(1) 
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Figure 2. Sample fits of the stellar kinematics of the nuclear region of NGC4051 using pPXF. Top left: fit of the nuclear spectrum; top 
right: fit of the spectrum at 1'.'5 E of the nucleus; bottom left: fit of the spectrum at the location corresponding to the blue turnover of 
the rotation curve at 1'.'2 SE of the nucleus; bottom-right: fit of the spectrum at the position of the red turnover of the rotation curve 
at 1'.'2 NW of the nucleus. The observed spectra are shown in black, the fits in red, the residuals in green, while in blue we show the 
spectral range not included in the fit of the central region. 



where a is a scale length, r is the radial distance in the plane 
of the galaxy, M is the mass inside r and G is the Newton's 
gravitational constant. Defining the coordinates of the kine- 
matical center of the system as {Xq, Yq), the observed radial 
velocity at position {R, *!/), where R is the projected radial 
distance from the nucleus in the plane of th e sky and ^ is 
the c orresponding position angle, is given by l|Barbosa et al.l 
I2OO6I ): 



Vr = Vs + 



R'^GM 



sin(i)cos(* - ^q) 



+ ^2)3/2 



COs2(>I> - ^o) + 



cos-^ (i) 



3/4 



(2) 



where Vs is the systemic velocity, i is the inclination of the 
disk (i = for face on disk) and \['o is the position an- 
gle of the line of nodes. The relations between r and R, 
and between a and A are: r = aR and a = aA, where 



a = Wcos2(vE' - ifo) + 



cos(i) 



The equation above con- 



tains six free parameters, including the kinematical center, 
which can be determined by fitting the model to the obser- 
vations. This was done using a Levenberg-Marquardt least- 
squares fitting algorithm, in which initial guesses are given 



for the free parameters. As the inclination of the disk is 
tightly coupled with M as Vj? oc Msin(z), it cannot be left 
as a free parameter. We have adopted the value of i = 41.4°, 
an estimate obtained from cos(i) = ^, where b and a are the 
semi-minor and semi-major axis of the large scale disk as 
quoted in the NASA/IPAC Extragalactic Database (NED) 
for this galaxy. 

The parameters derived from the fit are: the systemic 
velocity corrected by the observatory motion relative to the 
local standard of rest Vs =716± llkms"\ *o =120° ± 1°, 
M =7.7±0.6 X lO^M© and A =39.7 ± 2.7 pc. The derived 
kinematical center is very close to the peak of the continuum 
emission, with Xo=4.9±1.4pc and yb=3.1±1.2 pc, where Xq 
and Yo are measured in relation to the location correspond- 
ing to the peak of the continuum. 

In Fig. [5] we present the derived rotation model in the 
left panel and the residuals of the stellar velocity field (ob- 
served minus modelled) in the right panel. We conclude that 
the stellar velocity field is well described by the Plummer 
potential - the residuals are close to zero over most of the 
IFU field. The highest residuals (~ 20kms"^) are observed 
within ~0'.'3 from the nucleus, where the stellar kinemat- 
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Figure 3. Stellar kinematic maps obtained from tlie pPXF fit. Top; radial velocity (left) and velocity dispersion (right) maps. Bottom: 
/i3 and /i4 Gauss-Hermite moments. The mean uncertainties are 10 kms~^ for radial velocity, 8 kms^^ for cr, and 0.03 for and h^. 
The dashed lines show the position of the line of nodes. 
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Figure 5. Rotating disk model for the stellar kinematics of NGC 4051 (left) and residual map - observed velocity field less model (right). 
The dashed lines marks the position of the line of nodes. 



ics fitting may have been affected by emission by dust and 
emission lines. 



4 EMISSION-LINE FLUX DISTRIBUTIONS 

We have fitted Gaussians to the observed emission-line pro- 
files in order to obtain the integrated fiux, radial velocity 
(from the central wavelength of the line) and velocity dis- 
persion (from the width of the line) . The corresponding flux 
maps are shown in Fig.\6\ The H2 A =2.0338, 2.1218 and 
2.2235 fim flux distributions are presented in the top left, 
middle and right panels, respectively, with mean uncertain- 
ties of 16%, for the flrst one, 5% for the second one and 9% 
for the third. Black regions identify locations where the line 
fitting failed due to low signal to noise ratios. The molecu- 
lar hydrogen emission is extended over most of the observed 
field. The highest flux values are observed at the nucleus, de- 
fined as the location of the peak of the continuum emission. 
The H2 distribution is extended towards the NE, between 
the direction of the radio axis (adopted as the one which 
connects the two radio peaks) and the Une of nodes. The H2 
fiux distribution shows als o a good agreement with t he [O ill] 
narrow band image from [Schmitt fc KinnevI (|l996l ). whose 
contours are overplotted in green on the top right panel of 
Fig.El 

The flux map in the Br7 narrow component is shown 
in the bottom left panel of Fig.|6] In the central region 
(r < 25 pc) we have fitted two Gaussian components in or- 
der to separate the narrow and broad emission line contribu- 
tions, while in regions further away from the nucleus a single 
Gaussian was enough for the fit. The signal-to-noise ratio in 



the top ~ 1" of the observed field was not high enough to 
allow the measurement of this line and this is why the region 
appears black in Fig.|6] The mean fiux uncertainty for Br7 
is 24%. The Br7 flux distribution has the highest fiux values 
at the nucleus and is approximately symmetric around the 
nucleus, not showing any elongation as observed in the H2 
emitting gas. 

The flux map of the [Cavill] A2.3211 ^m coronal emis- 
sion line is presented in the bottom middle panel, where 
the mean flux uncertainties are 15%. The [Cavill] emission 
peaks at the nucleus and is resolved, extending up to O'.'S 
from the nucleus, which corresponds to a projected distance 
of 36 pc at the galaxy. 

The line ratio H2 A2.1218 /im/Br7 (only narrao compo- 
nent for Br7) can be used as a diagnostic for the excita- 
tion mechanism of th e molecular hydrogen emission lines 
(eg. Riffel et al.l |2006| ). with higher H2/Br7 ratios being 
interpreted as a larger contribution from shocks or X-ray 
to the H2 excitation. We present this line ratio map in 
the bottom-right panel of Fig.[5] The lowest values are 
H2 A2.1218 ^m/Br7 «1 observed at the nucleus and to the 
N. The highest values reach H2 A2.1218 /im/Br7 ~8 and are 
observed predominantly in two regions, one aproximately 
I'.'O W of the nucleus and another at O'.'S E of the nucleus. 
We note that these two regions are close to the tips of the 
compact 3.6 cm radio structure (black contours in the top 
middle panel). Nevertheless, such a direct interpretation of 
this map should be considered with caution, as the flux dis- 
tributions and kinematics (see next section) are obviously 
different for the H2 and Br7 emission lines, implying that 
may originate in different regions of the galaxy. 
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Figure 6. Top) From left to right: H2 A=2.0335, 2.1218 and 2.2235 /im flux maps; Bottom) left: Br7 flux distributions; middle: 
[Ca Vin]A2.3211 /jm flux distributions and right: H2 A2.1218/Br7 line ratio. The thick black contours overlaid to the H2A2.1218^m 
intensity map are from the VLA radio 3.6 cm continuum image, thin black lines are isointensities contours for each panel and the green 
contours are from an [Olll] narrow-band image from HST. The spatial scale and orientation shown at the top-left panel are the same for 
all panels. The dashed white line represent the line of nodes of the stellar velocity field and the full white line represent the PA which 
connects the two radio emission peaks. 



5 GAS KINEMATICS 

In the left panels of Fig.[7]we present the radial velocity field 
obtained from the central wavelengths of the H2 A2.1218 yLxn 
and Br7 emission lines, with mean uncertainties of 4kms^^ 
and 9kms~^, respectively. We chose the H2 A2.1218 pim line 
to represent the H2 velocity as it is stronger and thus present 
smaller uncertainties in the measurements than the other H2 
emission lines. The systemic velocity of the galaxy, derived 
from the stellar kinematics modelling, has been subtracted 
from all the emission line velocity plots. 

The H2 velocity field shows a "rotation pattern" similar 
to that of the stars, with the NW side receeding and the SE 
side approaching, although it is quite clear that there are 



other important kinematic components, evidenced by large 
deviations from simple rotation. Particularly conspicuous is 
a blueshifted region to the NE, showing velocity values up 
to ~ — lOOkms"^, extending by ~l" from the nucleus. The 
Br7 velocity field shows no rotation and the total velocity 
range is only «50kms^^. 

In the right panels of Fig. [7] we present the velocity 
dispersion (ct) maps obtained from measurements of the 

FWHM of the emission lines, such that cr = E^^m. The a 

' 2.355 

values were corrected for the instrumental broadening and 
the mean uncertainties are 6% for H2 and 22% for Br7. The 
H2 (T map has values in the range ~40-100 kms~^. A par- 
tial ring of low velocity dispersion values [a ~ 45kms~^) 
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is observed surrounding the nucleus, while higher values 
(80 ^ a < lOOkms"^) are observed over most of the re- 
maining field. The Br7 a map presents the highest values 
of up to « lOOkms"^ at the nucleus - but there are un- 
certainties in these a values due to the broad component 
contribution - and at «0.8 arcsec to the SW, close to the 
western tip of the 3.6 cm radio structure. At positions away 
from the nucleus the Br7 a is lower than 40kms~^. 

5.1 Gas "Tomography" 

The high spectral resolution of the data has allowed us 
to slice the emission-line profiles of H2 A2.1218 fim and 
[Ca vill] A2.3211 /xm into a sequence of velocity bins. With 
this "tomography" technique, we can sample the kinematics 
along the whole emission- line profile, including the wings. 
In order to obtain the "tomography" images, we resampled 
the spectra into bins of 1 A with the scombine iraf task 
and then combined the two sets of observations into a single 
data cube using the tasks SCOMBINE and imcombine. The 
velocity slices were obtained after subtraction of the contin- 
uum determined as averages of the fluxes from both sides of 
the emission line. The slices correspond to velocity bins of 
« 42kms"^ (3 A) and are shown in Figs. [8] and [O] In these 
figures, each panel presents flux levels in logarithmic units 
for the velocity slice shown. The zero velocity is adopted as 
the systemic velocity of the stars obtained from the stellar 
kinematic modelling. 

The slices trace the gas from negative (blueshifts) to 
positive (redshifts) velocities relative to the systemic veloc- 
ity of the stars. For H2, the highest blueshifts, which reach 
« — 190 kms~^, are observed predominantly to the NE along 
a curved elongated structure similar to the one observed 
in the radial velocity map (top-left panel of Fig. [7]). This 
structure - whose morphology in Fig. [S] can be described as 
a curved arm - dominates the emission from « —190 to 
~ — 60kms~^. At velocities close to systemic, the emission 
is approximately symmetric and dominated by a region of 
~ 1" radius centered at the nucleus. As the slices reach pos- 
itive velocities, the dominant structure is another curved 
arm extending to «2" to the W of the nucleus. The highest 
redshifts reach ~ 160kms~^. 

For [Cavill], we show in Fig.[9]the highest velocity bins 
and exclude a few low velocity bins where the kinematics are 
similar. The highest blueshifts reach ~-250kms~^, which 
are higher than those observed for the H2 emitting gas, while 
the highest redshifts reach velocities of 170 km s^^ similar to 
the ones observed for H2. 

We opted not to show slices in the Br7 emission-line 
profile due to the fact that its narrow component is indeed 
narrow and quite symmetric, and close to the nucleus it is 
hard to deblend it from the broad line profile, which intro- 
duces too much uncertainty in the derived kinematics. 



6 DISCUSSION 

6.1 Stellar Kinematics 

As observed in the top-left panel of Fig.Oand in Fig.Blthe 
turnover of the rotation curve occurs at only «55 pc from the 
nucleus, suggesting that the stellar motions are dominated 



by a highly concentrated gravitational potential. This is also 
supported by the small value obtained for the scale length 
{A — 39.7 pc) from the modelling of the velocity field. The 
derived paramet ers are in appr o ximat e agreement with the 
ones derived by iBarbosa et al.l (I2OO6I ) from a similar mod- 
elling using optical IFU data for a field-of-view of ~ 7" x 5" 
(315x225 pc^). The exception is our value of = 120°, 
which is ~ 13° larger than theirs. On th e other hand, our 
is ~ 13° smaller than the one derived bv lDumas et al.l (|2007l ) 
from stellar kinematics obtained over the much larger field- 
of-view of « 50"x40" (2250xl800pc^). We atribute these 
differences to the differences in fleld-of-view, considering also 
that our field-of-view, which corresponds to 130xl80pc^ at 
the galaxy, is only sampling the stellar kinematics very close 
to the nucleus. 

For an orientation of the line of nodes of 'I'o — 120°, 
and assuming that the spiral arms of NGC 4051 are trailing 
(orientation of the arms is shown in Fig.[T|, we conclude that 
the NE is the far side and the SW is the near side of the 
galaxy. 

The stellar velocity dispersion values show patches with 
lower values (~40-50kms~^) on top of a common back- 
ground of values ranging from 60 to 70kms~^. A possi- 
ble interpretation is that these patches are colder regions 
with more recent star formation than the underlying bulge. 
This is supported by optical spectra of the nuclear re- 
gion which show clear sign atures of intermediate age stars 
(jCid Fernandes et al.|[200a) . A very small a drop is observed 
right at the nucleus, but we take this result with caution 
because of the obvious contamination of the spectra of the 
nuclear region by emission from dust and broad emission 
lines. 

We can estimate the mass of the SMBH (Mbh) from the 
bulge stellar velocity dispersion (a*) as \og{MBH /Mq) = 
a + /31og((7./(jo), wh ere a = 8.13 ± 0.06, /3 = 4.02 ± 0.32 
and (To — 200kms~^ (|Tremaine et al. I I2OO2D . Adopting cr. ^ 
60kms^^ as representative of the bulge (top right panel 
of Fig.EI, we obtain Mbh = 1.1 ± 0.3 x IO'^Mq. This 
value is in good agreement with those obtained by previous 
authors from reverberation mapping a nd scaling relations 
(|Shemmer et al.ll2003l : fKaspi et al.ll2000l ). For this mass, the 
radius of influence of the SMBH is ~ 1.3 pc, thus not re- 
solved at the spatial resolution of our data. 

6.2 Gas Kinematics 

The simultaneous observation of the stellar and gaseous 
kinematics allowed us to construct residual maps for the 
gaseous kinematics relative to the stellar kinematics model 
described in section l3^ The residual map for the H2 emit- 
ting gas is presented in Fig. llUI The most conspicous feature 
in this map is the elongated structure to the NE, which 
shows blueshifts of up to ~ — lOOkms"^, also seen in the 
corresponding velocity slices in Fig.|Sl These blueshifts could 
be either due to an inflow, if the gas is located in the plane, 
as the NE is the far side of the galaxy; or to an outflow, if 
the gas is extended to high latitudes (above the plane), e.g. 
in a conical structure oriented towards us. We discuss below 
these two possibilities. 

The outflow interpretation is supported by the follow- 
ing facts: (1) that the elongated NE structure is approxi- 
mately oriented along the radio jet; (2) a similar structure 
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Figure 7. Left: Velocity field for the H2 A2.1218 /xm (top) and Br7 (bottom) emission lines. Right: Velocity dispersion maps for the same 
emission lines. The contours and lines are as described in Fig.|6] as well as the spatial scale. 
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Figure 8. Velocity slices along the H2 profile with a velocity bin of ^ 42kms~^. The cross marks the position of the nucleus and the 
intensity scale is logarithmic. The zero velocity is the systemic velocity obtained from the stellar kinematics modelling. The near and far 
side of the galaxy are indicated in the top left panel (under the assumption that the spiral arms are trailing). 



12 Rogemar A. Rijjel et al. 



CM 



O 
<D 
CO 
O 



CM 




1 -1 

(arcsec) 

Figure 9. Same as Fig. |8] for the [Cavill] emission-line profile. 



and blueshifts have been observed bv iBarbosa et all l|2008l . 
-hereafter BOS) in IFU observations of [S ill] A9069 emitting 
gas, who found, in addition, redshifts to the SW, which ap- 
pear to be due to a counterpart conical outflow related to 
the SW part of the radio-jet, probably located behind the 
galactic plane and being directed away from us; (3) the ve- 
locity dispersion and line ratio maps which show increased 
values in regions close to the tips of the radio jet. 

The inflow interpretation is, on the other hand, sup- 
ported by: (1) the fact that the orientation of the NE 
blueshifted region is not well aligned with the direction of the 
radio jet, but is shifted by ~0'.'5 to the S-SE (see Fig.|S}; (2) 
the fact that the blueshifted region is elongated and curved 
as if it belonged to a spiral arm ending at the nucleus, while 
the blueshifted structure observed in the [S ill] A9069 emit- 
ting gas by BOS shows a more conical shape whose axis is 
better aligned with the direction of the radio jet; (3) the fact 
that there is no redshifted counterpart in H2 as observed for 
the [S III] emitting gas. There is instead a redshifted struc- 
ture observed in the velocity slices (see Fig.|8]) also curved 



as belonging to a spiral arm in the near side (SW) of the 
galaxy. 

Considering the arguments above we favour the inflow 
interpretation. The blueshifts to the NE and redshifts to 
the SW observed in the H2 emitting gas would then be due 
to inflow of molecular gas in the galaxy plane along nu- 
clear spiral arms. The presence of such arms is supported 
by the struc ture map of the nuclear r egion of NGC 4051 
presented bv lSimoes Lopes et al.l (|2007| ). which shows simi- 
lar spiral structure with the same curvature and orientation 
of the arms in an HST optical image of the nuclear region of 
NGC 4051. The structure map also shows a lot of obscura- 
tion against the near side of the galaxy, which would explain 
why the flux of the redshifted gas emission is fainter than 
that of the blueshifted gas on the far side of the galaxy. The 
enhancements in the gas velocity dispersions and increase 
in line emission ratios observed in H2 and Br7 emitting gas 
can be understood if the radio jet is launched at a relatively 
small angle to the galaxy plane compressing the circumnu- 
clear ISM close to the nucleus. 
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In a recent study aimed at investigating the interac- 
tion between the radio jet and line emission in the NLR 
of the Seyfert galaxy ES0 428-G14 using GNIRS IFU data 
we found a close association of the [Fell] and Hi emission 
line distribu tions and kinema tics with the radio emission 
distribution dRiffel et al.ll2006l). a result also found in some 
studies (eg. Falcke. Wilson fc SimpsonI 1 19981 : iBicknell et al.l 



(eg; 



I2OOOI: iTilak et al 



Kaiser et al 



2005 , BOS), but also c ontested by others 



2000l : iDas eral]l2005l ). In ES0 428-G14, 



iRiffel et al. (2006) concluded that the H2 kinematics was 
distinct from that in the other emission-lines: although part 
of the II2 emitting gas was affected by interaction with the 
radio jet, most of it was located in the galaxy plane and was 
not affected by the radio jet. A less disturbed kinematics for 
H2 than for other emission lines was also found in a near-IR 
study of the gas emission around the nuclei of the Seyfert 
galax ies NGC 2110 and Circicus (jStorchi-Bergmann et al.1 
I1999I ). These studies thus also support our favored inter- 
pretation in the present study: most of the H2 emission is 
located in the plane and is not related to the outflow. It is 
instead, flowing towards the nucleus along spiral arms. 

The Br7 emitting gas is almost all at the systemic veloc- 
ity as observed in the bottom-left panel of Fig. [7] and seems 
to show no rotation, which leads us to conclude that most 
of the Br7 emitting gas is not restricted to the plane of the 
galaxy. 



6.3 [Cavill] coronal emission 

Coronal lines are forbidden transitions from highly ionized 
species which extend from the unresolved nucleus up to dis- 
tances between a few tens to a few hundreds of parsecs, 
and usually present blue wings and are broader than low 
ionization lines (e.g. iRodrfguez-Ardila et al.ll2006l . and ref- 
erences therein). Thus, it is possible to resolve the coro- 
nal line region for the closest active galaxies. In NGC 4051, 
the [Cavill] coronal emission line at 2.3211 /xm is indeed 
resolved, presenting extended emission in a circular region 
around the nucleus with diameter of «75pc. Blueshifts of 
up to — 250kms~^ are observed very close to the nucleus 
(within the inner «25pc), a value that is much higher than 
that observed for Br7 and II2 emission lines. The mean ve- 
locity dispersion is 150kms~^, which is about 2-3 times 
higher than that obtained the other lines as well. This kine- 
matics together with the compact flux distribution supports 
the interpretation that the coronal lines are produced close 
to the nucleus, and probably in the transition region between 
the BLR (broad-line region) and the NLR as suggested by 
previous authors. The spatial extent of the emission, the 
blueshifts and the a values observed in NGC 4051 are sim- 
ilar to the ones obtained for the Seyfert 2 galaxy Circinus, 
while for more active galaxies such values tend to be higher 
l|Rodn'guez-Ardila et al1l2006l ). 



6.4 Mass of the emitting gas 

The mass of hot II2 can be estimated as (e.g. IScoville et al.l 
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Figure 10. Difference between the H2 and stellar kinematics. 
The white dashed line shows the orientation of the line of nodes 
and the full white line shows the PA of the radio emission. 



5.0776 X 10' 



F, 



ff2-'*2.1218 



erg s 



cm" 



Mpc 



(3) 



where rup is the proton mass, i<if2A2.i2i8 is the line 
flux, d is the galaxy distanc e and Mh^, is give n in 
solar masses. As proposed by IScoville et all (|1982| ). we 
assume a typical vibrational temperature of Tvib = 
2000 K, which implies a population fraction f^^i j^^ = 



1.22 X 10" 



and 



a transition probability ^.s ni = 

(|Turner. Kirbv-Docken fc Dalgarnclll977l ). 
The line flux should be corrected for the intrinsic 
E(B - V), but for NGC4051, E(B - V) is smaU 
(|Riffel. Rodriguez- Ardila fc Pastorkall2006l ) and thus we do 
not apply any correction. For the E-NE inflow region we 
obtain -f£r2A2.i2i8 = 2.1 x 10~^^ ergs~^ cm~^ integrated 
within the approximately triangular region which includes 
blueshifts larger than — lOOkms"^ in Fig. 1 101 and obtain a 
mass of Mh2 ~ 9.3 Mq. Integrating now the H2 A2.1218 
emission over the same fleld where the Br7 is observed 
(the central 130xl35pc^) we obtain Fjj2A2.1218 = 1-5 x 
10"" erg s"^ cm" ^ and thus ^ 66M0. 

The above masses are small, but we point out that in 
the nuclear region of galaxies it has been shown t hat the hot- 
to-co ld mass ratio ranges between lO"'^ to 10~^ l|Dale et al.l 
l2005h . suggesting that the total mass of molecular gas can be 
orders of magnitude larger than the value measured directly 
from the H2 emission. 

The mass of the ionized hydrogen can be estimated 
as Mhii ~ rUp NeVnii, where A^e is the electron density 
and Vhii is the volume of the e mitting region. Replacing 
A^e Vhii by the expression given in IScoville et al.l l|l982l ) we 
obtain: 
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where Mhii is given is solar masses, and we assume an 
electron temperature T — 10* K and electron density A^'e = 
100 cm~^. The total integrated Br7 flux is Fsr-y ~ 5.8 x 
lO'^'^ergs-^cm"^ and thus Mhii « 1.4 x lO'^M©. The 
mass of ionized gas is thus about 2000 larger than the mass 
of hot H2 gas. 

Using the above calculated mass for the E-NE 

blueshifted region we can estimate the flux of inflowing hot 
molecular gas under the assumption that this gas is stream- 
ing towards the nucleus. We calculate the hot H2 mass 
crossing a circular cross section by Mh2 = 2mp Nh2 virr^ , 
where v is the inflowing velocity and r is the radius of 
the circular cross section. The maximum blueshift veloc- 
ities occur at l'.'2 E-NE of the nucleus where the cross- 
section radius of the structure shown in Fig. [10] is esti- 
mated to be r =0'.'55« 25 pc. Assuming v = lOOkms"^ 
and conical geometry to calculate the H2 density, we obtain 
= 2.5 X lO^^gs^^ Ki Ax lO"^M0yr-^ This rate is 
probably up to 2 times larger if we consider a similar inflow 
along the W-SW spiral arm, observed in redshift. 

The mass accretion rate necessary to power the active 
nucleus can be estimeted by M = where Lboi is the 

bolometric luminosity, c is the light speed and r] is the effi- 
ciency of conversion of the rest mass energy of the accreted 
material into radiation. Assuming a typic al value of ?? » 0.1 
and using Lboi = 2.7 x lO"^ ergs"^ from lOgle eVsI\ ()2004l l 
we obtain M ^ 4.7 x 10~^ Mq/jt. The H2 mass inflow rate 
calculated above is then about 100 times smaller than the 
accretion rate needed to power the AGN of NGC4051, sup- 
porting additional contribution of cold non-emitting gas to 
the mass inflow rate. 



6.5 H2 excitation 

The H2 lines can be excited by three distinct processes: 
fluorescent excitation through absorption of soft-UV 
photons (912-1108 A) in th e Lyman and Werner bands 



Black fc van Dishoecklll98'it). e xcitati on by X-ray heating 



.,Malonev . HoUenbach fc TielensI Il996l ) and excitation by 
shocks (Hollcmbacli & McKed Il989l ). The first is usu- 
ally considered a non-thermal process while the other 
two are considered thermal processes. Several stud- 
ies investigated th e H2 e xcitation mechanisms ( eg. 



1997 


: iQuillen et al. Il999l: 


Bellamy fc Tadhuntcr 


2004 


: iRodriEuez-Ardila. Riffel & Pastoriza 


I2OO5I: 


Riffel. Rodriguez- Ardila & Pastorizal 20061: 


Riffel et all 



IQuillen et al.l l|l999i) have looked for correlations of 
the H2 emission with radio 6 cm and hard X-ray fluxes. 
They found no correlation with X-rays, suggesting X-rays 
heating is not the dominant H2 excitation mechanism, 
and found a weak correlation with radio 6 cm, suggest- 
ing that no single mechanism is likely to be responsible 
for the molecular hydr ogen excitation in Seyfert galaxies. 
iDraine fc Woodd i 1990|) propose that most of the H2 line 
emission originates in molecular gas which is heated by 
transient X-ray irradiation. H ard X-rays fro m the AGN 
also have been proposed by IWilman et. all (|200G| ) and 



[Bellamy fc Tadhuntej (|2004l ') as the dominant excitation 
mechanis m of H2 emission in Cygnus A. Zuther et al.l (|2007f ) 
usmg thelMa lonev. HoUenbach fc TielensI l| 19961 ) X-rays ex- 
citation models found that the X-ray emission can account 
for some H2 excitation for the starburst/Seyfert galaxy 
Mrk609. For the case of ES0428-G14 we concluded that 
shocks produced by a radio jet play a fundamental role in 
the kinematics and excitation of the H2A2.I2I8 emission 
line in regions co-spatial with the radio emission and X- 
rays are importan t in regions away from the radio structures 
iRiffel et al.ll2006l) . 

The H2 A2.2477/A2.1218 line ratio is commonly used 
to distinguish between thermal (~0. 1-0.2) and fluores- 
cent (~0.5 5) exci t ation mechanisms for the H2 emit- 
ting gas |Mouril 1994 : Reunanen. Kotilainen fc Prietol 



l2002l : [Rodriguez- Ardila. Riffel fc Pastorizall2005l ). while the 
H2 A2.0338/A2.2235 ratio is used to determine the thermal 
excitation temperature. The H2 A2.2477 /im emission line is 
only marginally detected in our individual spectra, so we 
integrated the flux of this line over the whole IFU field ob- 
taining H2 A2.2477/A2.1218 « 0.12±0.02 indicating that the 
emitted H2 is excited by thermal processes. 

The rotational and vibrational temperatures for the 
molecular gas can also be used to distinguish between ther- 
mal and fluorescent excitation: for thermal excitation the 
vibrational and rotational temperatures must be similar and 
for fluorescent excitation the vibrational temperature must 
be higher than the rotational temperature - non-local UV 
photons overpopulate the highest energy levels compared 
to the population distribution expected for a Maxwell- 
Boltzmann distribution. The rotational temperature can be 
obtained by T^ot = -1113/ln(0.323 ^«2A2.o338 -^ ^j^^ 

^ H '2 A 2 . 2235 

brational temperature by Tvib = 5600/lnfl.355 £|£2i2:i2ii') 

t ' ^ ^jr2A2.2477 ' 

(jReunanen. Kotilainen fc Prietol I2OO2I ). For NGC4051 
we obtain H2 A2.0338/A2.2235 = 1.8 ± 0.28 and 
H2 A2.1218/A2.2477 = 8.4 ± 1.02, thus T^ot = 20521^1 K 
and Tvib = 2305^120 "^^^ similarity of these two values 
reinforce thermal processes as the dominant excitation 
mechanism of the molecular hydrogen. The value obtained 
for the vibrational temperature supports the assumption 
(Tvib = 2000 K) used in equation[3| to obtain the hot H2 
mass. 

In order to test if X-ray emission can account for the 
excitation of the H2 A 2. 1218 /xm emission line for NGC4051 
we ha ve used the models of iMalonev. HoUenbach fc TielensI 
l|l996t ) to estimate the emergent H2 flux of a gas cloud illu- 
minated by a source of hard X-rays with an intrinsic lumi- 
nosity Lx- Th e cooling is given by the effective ionization 
parameter ^efi ([Zuther e t al. 2007j): 



^off = 1.26 X 10" 



fx 



(5) 



where fx is the incident hard X-ray flux at the distance 
d[pc] from the X-ray source, n5[10~*cm~^] is the total 
hydrogen gas density, and Af22 [lO^'^cm^] is the attenuat- 
ing column density. Using /x=2.3x 10~^^ ergcm"^ s~^ 
and iVi/=1.32xl0^° cm~^ obtained from the 
ASCA Tartarus database w e ca n calculate ^cff- 
[Malonev. HoUenbach fc" TielensI { 19961 ) calculate emer- 
gent fluxes for two gas densities, lO^cm"^ and lO'^cm""^. 
We calculate the effective ionization parameter for three 
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Table 1. Comparasion of the observed H2 A2.1218 f im fluxes and 
calculated using models of Malone v. HoUenbach fc T iclcns (1996) 
for an aperture of O'.'lxO'.'l for hydrogen densities n = 10^ cm""^ 
and n = 10'^ cm"'' . 





Observed 


n = 


10^ cm ^ 


n = 10^ 


cm ^ 


d 


logCFffJ 


log(Ceff 


) log(i^Hj 


log(€eff) 


logCi^ilJ 


25 


-16.4 


-1.7 


-14.7 






50 


-17.0 


-2.3 


-17.2 


-0.3 


-15.8 


75 


-17.3 


-2.7 


-17.0 


-0.7 


-15.7 



different distances of the A GN (25, 50 and 75 pc), and for 
the sa me gas densities of iMalonev. Hollenbach fc Tielend 
l|l996t ). Using this and their Fig. 6(a,b) we can obtain the 
emergent intensity of the H2 A2.1218 /xm emission line. 
In Table[l] we present the calculated effective ionization 
parameter and the emergent flux for an aperture of 
O'.'lxO'.'l - corresponding to a solid angle 2.34x10"^^ sr - 
in ergcm~^s~^. In this table we also present the observed 
H2 A2.1218 ^m obtained by the average of the H2 flux in a 
ring with radius d and width of 0'.'2. 

By comparing the emergent fluxes from Table[T] ob- 
tained by X-ray heating models with the observed fluxes 
we conclude that excitation by X-ray heating can account 
for most of the observed H2 flux. 

Another line ratio commonly used to investigate the 
H2 excitation is H2 A2.1218/Br7. In starburst galaxies, 
whe re the main heating agent is UV rad i ation , H2/Br7 < 
0.6 iRodn'guez-Ardila. Riffel fc Pastorizal I2005D . while for 
Seyferts this ratio is larger (0.6 <H2/Br7 < 2) because 
of additional H 2 emission excited by shocks or by X-ray s 
from the AGN (|Rodrfguez-Ardila. Riffel fc Pastorizall2005l ). 
As observed in the bottom-right panel of Fig.|6l NGC 4051 
presents H2/Br7 = 1 ± 0.2 at the nucleus, which is a typical 
value for AGNs. The highest values reach H2/Br7 = 8 ± 1.9 
and are observed E-NE and to the W-SW of the nucleus, 
close to the tips of the radio jet, supporting some contri- 
bution of shocks by the radio jet to the excitation of the 
molecular hydrogen in these regions. On the other hand, the 
interpretation of this line ratio should be considered with 
caution due to the fact that the bulk of II2 emission seems 
to originate from the disk while the Br7 emission seems not 
to be restricted to the disk. 



7 CONCLUSIONS 

We have analysed two-dimensional near-IR if— band spec- 
tra from the inner « 150 pc of the Narrow Line Seyfert 1 
galaxy NGC 4051 obtained with the Gemini NIFS integral 
field spectrograph at a sampling of 4.5x4.5 pc^ at the galaxy 
and spectral resolution of ~3 A. We have mapped the stellar 
and gaseous kinematics, and the emission-line flux distribu- 
tions and ratios of the molecular and ionized hydrogen. The 
main conclusions of this work are: 

• The turnover of the stellar rotation curve is at only 
~55pc from the nucleus, suggesting that the stellar mo- 
tions are dominated by a highly concentrated gravitational 
potential, a result confirmed by modeling using a Plummer 
gravitational potential, as we obtain a sma ll scale length, 
A «39 pc. This result supports the findings of lBarbosa et al.l 



(|2006l ) based of optical data. The mean velocity dispersion of 
the bulge (~60kms"^) implies a SMBH mass of ~ 10'^ M©. 
Within the bulge, we find patches of lower velocity disper- 
sions, which we attribute to recent star formation. 

• The gas kinematics is distinct for each emission line: 
the Br7 emission line shows velocities restricted to within 
~30kms~^ from systemic, without evidence of bultk mo- 
tions, suggesting that the ionized gas is not restricted to the 
galactic plane, while much larger blueshifts and redshifts 
are observed for both the [Ca vill] and H2 emitting gas. The 
Br7 velocity dispersion is small (~40kms~^) over most of 
the field, except around the nucleus at a location close to the 
SW tip of the radio jet, suggesting some energy injection in 
the gas by interaction with the radio jet. 

• The [Cavill] coronal emission line is compact but re- 
solved, extending over a circular region of radius ~35pc 
around the nucleus. It also present the largest ve- 
locity blueshifts (— 250kms~^) and velocity dispersion 
(150 kms"^) of the emission lines studied here, supporting 
an origin close to the active nucleus, possibly in the transi- 
tion region between the BLR and the NLR. 

• The H2 emitting gas seems to be mostly restricted 
to the galaxy plane. The most conspicuous kinematic fea- 
tures are a curved elongated blueshifted structure to the 
NE, interpreted as gas infiow along a nuclear spiral arm 
in the far side of the galaxy, and a curved, redshifted 
structure to the SW, interpreted as gas inflow along a 
nuclear spiral arm in the near side of the galaxy. Esti- 
mates of the mass inflow rate in hot H2 gives ~ 
4 x 10"^MQyr"\ which is ~ 100 times smaller than the 
nuclear accretion rate necessary to power the active nu- 
cleus of NGC 4051, supporting the presence of additional 
inflow of cold non-emitting gas. This is not the first time 
we have been able to map infiows towards an active nu- 
cleu s along nuclear spiral arms. In two prev i ous s tudies 
we (|Fathi et al.l I2OO6I : IStorchi-Bergmann et al] bOOTl ) have 
mapped streaming motions towards the active nucleus along 
nuclear spirals in the galaxies NGC 1097 and NGC 6951 
although in ionized gas. As nuclear spirals are ubiqui- 
tous around active galactic nucl ei (Sin ioes Lopes et al ] |2007l : 
IPrieto. Macieiewski fc ReunanenlbOoH ). such infiows seem 
to be an "universal" mechanism to bring gas inwards to 
feed the SMBH within the inner few hundred parsecs of the 
galaxy. 

• The total mass of hot H2 is estimated to be of the order 
of 66 M©, while that of Hll is estimated to be 1.4 x lO'^ Mq. 

• From the H2 line ratios we conclude that H2 is excited 
by thermal processes - heating by X-rays from the AGN and 
shocks produced by the ra dio jet. We conclude, based in 
X-ray excitation models of iMalonev. Hollenbach fc Tieleng 
(| 19961 ) that X-ray heating can account for the observed emis- 
sion, but the H2 A2.1218 /im/Br7 line ratio supports some 
contribution from shocks in the regions where the radio jet 
interacts with the H2 emitting gas. 
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